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Abstract 

The development of accurate theoretical failure, fatigue, and wear models for ultra-high molecular weight polyethylene 
(UHMWPE) is an important step towards better understanding the micromechanisms of the surface damage that occur in load 
bearing orthopaedic components and improving the lifetime of joint arthoplasties. Previous attempts to analytically predict the clin- 
ically observed damage, wear, and fatigue failure modes have met with limited success due to the complicated interaction between 
microstructural deformations and continuum level stresses. In this work, we examined monotonic uniaxial and multiaxial loading to 
failure of UHMWPE using eight failure criteria (maximum principal stress, Mises stress, Tresca stress, hydrostatic stress, Coulomb 
stress, maximum principal strain, Mises strain, and chain stretch). The quality of the predictions of the different models was assessed 
by comparing uniaxial tension and small punch test data at different rates with the failure model predictions. The experimental data 
were obtained for two conventional (unirradiated and gamma radiation sterilized in nitrogen) and two highly crosslinked ( 1 50 kGy, 
remelted and annealed) UHMWPE materials. Of the different failures models examined, the chain stretch failure model was found 
to capture uniaxial and multiaxial failure data most accurately for all of the UHMWPE materials. In addition, the chain stretch 
failure criterion can readily be calculated for contemporary UHMWPE materials based on available uniaxial tension data. These 
results lay the foundation for future developments of damage and wear models capable of predicting multiaxial failure under cyclic 
loading conditions. 

© 2004 Orthopaedic Research Society. Published by Elsevier Ltd. All rights reserved. 
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Introduction 

Mechanical failure of ultra-high molecular weight 
polyethylene (UHMWPE) components is one of the 
major factors controlling the lifetime of total joint artho- 
plasties. The failures that occur in vivo are often intro- 
duced by contact damage in the form of wear and 
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fatigue [19] caused by relative sliding motion between 
different components, primarily due to metal-to-polymer 
interactions. To increase the long-term performance of 
total joint arthroplasties, investigators have attempted 
to understand the micromechanisms underlying surface 
damage in UHMWPE components [27,28]. The reason- 
ing for this approach is that a fundamental understand- 
ing of the micromechanisms of surface damage will 
facilitate predictions of performance on the component 
level, which in turn can enable design and material 
improvements. Despite these efforts, the micromechan- 
ics controlling the evolution of surface damage in 
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UHMWPE components are still not fully understood, 
and the state of the art in modeling the failure events is 
far from complete. 

Modeling UHMWPE damage and failure, based 
on micromechanical events, is a challenging task. On 
the local scale, the material consists of intermixed 
amorphous and semicrystalline domains. The material 
behavior is further complicated by the evolving micro- 
mechanical state due to applied external loading. Failure 
of UHMWPE on the microscale is driven by molecular 
chain scission [12,17]. The progressive failure events are 
promoted by the magnified stress and strain state in the 
neighborhood of already broken molecular chains, 
resulting in damage localization that is manifested by 
growth and coalescence of microvoids (also referred to 
as cavitation) [24,29]. 

In general, failure of any polymeric component can 
be defined as the inability to perform an intended func- 
tion. Examples of potential failure modes are excessive 
plastic deformation, static fracture, fatigue fracture, 
and wear. Our study focused on static fracture, defined 
as loss of material continuity through the formation of 
new surfaces within the UHMWPE component. This 
failure mode is the most fundamental, and any failure 
model of common clinical failure modes, such as fatigue 
or wear, will need to also address failure under mono- 
tonic loading conditions. 

There are two general classes of failure models for 
polymers: isostate models and evolving state models. 
An isostate failure model is a failure model in which 
the material state is assumed not to change during the 
applied deformation. The most general isostate failure 
model is given by 

/(T,E) >/m t , (1) 

where T is the current stress state (expressed, for exam- 
ple, in terms of the Cauchy stress tensor), E is the cur- 
rent strain state (expressed, for example, in terms of 
the logarithmic strain tensor), J{ ■ ) is a pre-specified fail- 
ure model, and / crit is a material parameter specifying 
the critical value of the failure condition. There are 
two main classes of evolving state models: state variable 
and history-based models. A general failure condition in 
state variable form can be written as 

/(£i>T,E) > fcriu (2) 

where & (j G [1, iV]) is a set of state variables specifying 
the current state. The state variables are assumed to 
evolve with the deformation according to the rate- 
equation: ti — g(£/,T,E), where g( • ) is a pre-specified 
evolution equation. 

An evolving state model in history-based form can be 
written as 

f /(T(/),E(/'))d/' >/ cri t, (3) 

J -OO 



where Eq. (3) symbolizes that the failure not only 
depends on the current values of stress and strain, but 
also on their complete time histories. One interpretation 
of the history-based integral model can be obtained by 
defining the integral on the left-hand side of Eq. (3) as 
the damage D{t), The rate of change of damage is then 
given by D(t) = /( T,E), and the history integral repre- 
sentation can consequently be seen as a special case of 
the state variable-based evolving state failure model. 

The goal of this study was to examine different 
approaches for predicting static fracture, and then in 
a later study use this necessary information to address 
fatigue and wear failures. In this work, we have used 
monotonic uniaxial and small punch (equibiaxial) tests 
to failure to examine the relative capabilities of different 
failure models to predict fracture for two conventional 
(unirradiated and gamma radiation sterilized in nitro- 
gen) and two highly crosslinked (150 kGy, remelted 
and annealed) UHMWPE materials. Currently, conven- 
tional UHMWPE is the dominant material for hip and 
knee total joint replacement implants. However, highly 
crosslinked UHMWPE materials are gaining in use, par- 
ticularly for hip replacement components, due to im- 
proved wear resistance [21]. 

Materials and methods 

A number of different types of models can potentially be used to 
predict failure of UHMWPE, including: stress-based models, strain- 
based models, energy-based models, linear elastic fracture mechanics 
(LEFM) models, and damage models [2]. The first four are constant 
state models, whereas the last type are evolving state models. An en- 
ergy-based model is unlikely to be predictive since large strain energy 
can be stored in hydrostatic compression but the material will not frac- 
ture in this mode. LEFM models are a specific case of a stress-based 
constant state model. Thus, in this study, we examined the following 
eight isostate failure models: maximum principal stress, Mises stress, 
Tresca stress, hydrostatic stress. Coulomb stress, maximum principal 
strain, Mises strain, and chain stretch. Most of these models are well 
known and are discussed in detail in most introductory solid mechan- 
ics textbooks [22]. For completeness, three of the lesser known models 
(the Coulomb model, the Mises strain model, and the chain stretch 
model) are discussed herein in more detail. 

The Coulomb stress is defined as 

^Coulomb = (ffl - (4) 

where cr, is the largest principal stress, <t 3 is the smallest principal 
stress, and p = 0.10 is taken as the friction coefficient. The Mises strain 
is defined analogously to the traditional Mises stress, i.e. 

^Mises = I( £l ~ fi 2)“ + ~ fi 3) 2 + - Cl) 2 ]* (5) 

where e, are the principal strains. 

The chain stretch failure condition is motivated by our recently 
developed physically-inspired constitutive model for UHWMPE, the 
Hybrid Model (HM) [9,10]. The HM is based, in part, on some of 
the constitutive modeling work by Boyce et al. [3,4,16] who success- 
fully modeled glassy amorphous polymers, semicrystalline polymers, 
and rubbers. However, these theories did not adequately describe 
UHMWPE; thus, the HM was developed and validated specifically 
for semicrystalline UHMWPE. A key feature of the HM is the concept 
of molecular chain stretch [5]. In the chain stretch failure model, an 
estimation of the average molecular chain stretch is given by the 
expression 
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where, as before, e 2 , and e 3 are the principal strains. This represen- 
tation enables a mapping from the applied deformation state to the 
local molecular stretch state and has been examined using molecular 
dynamics and Monte Carlo simulations [8]. In the chain stretch model, 
failure is assumed to occur when the average molecular chain stretch 
reaches a critical level. 

To evaluate the quality of the predictions of the different failure 
models, appropriate experimental test methods and well-characterized 
materials must be used. We examined four different types of 
UHMWPE obtained from ram-extruded rods of GUR 1050 that were 
subjected to different radiation and heat treatments. The materials [22] 
were: as-received unirradiated, gamma radiation sterilized in Nitrogen 
with 30 kGy, gamma radiation treated with 100 kGy (giving a highly 
crosslinked microstructure and subsequently annealed at 110°C 
for 2 h), and gamma radiation treated with 100 kGy (giving a highly 
crosslinked microstructure and subsequently remelted at 150 °C). 
The microstructures of these materials were examined elsewhere [22]; 
e.g., the degree of crystallinity of the three materials was determined 
to be 0.50 for the unirradiated, 0.51 for the sterilized material 
(30 kGy y-N 2 ), 0.61 for the crosslinked material heat treated at 
110°C, and 0.46 for the crosslinked material heat treated at 150 °C. 
The four materials were tested in uniaxial large strain monotonic ten- 
sion to failure and in small punch testing to failure. 




Fixture 



Fig. 1 . Finite element mesh used in the small punch simulations. In the 
testing, the spherical punch head indented the disc-shaped specimen 
that was held in place in the fixture by frictional forces. The figure 
shows the deformed mesh at a punch displacement of 2.0 mm. 



Uniaxial tension to failure experiments 

The stress-strain behavior from the uniaxial tension to failure tests 
was reported previously [16]. Uniaxial tension behavior of other types 
of conventional and crosslinked UHMWPE were reported in other 
studies [1,13,15,18,20,23,25]. Briefly, specimens were prepared and 
tested according to ASTM D638 [6]. Testing was performed using a 
4204 Instron electromechanical load frame, equipped with a 5 kN load 
cell, under displacement control at crosshead displacement rates of 30, 
75, and 150mm/min, approximately corresponding to initial true 
strain rates of 0.0 1/s, 0.02/s, and 0.05/s. Axial strain was measured 
using non-contacting video extensometer. A total of 60 specimens were 
tested (5 specimens per material x 4 materials x 3 rates). 

Small punch testing to failure 

The load-displacement behavior from the small punch failure tests 
was reported previously [19]. Briefly, the small punch testing was per- 
formed on disc shaped specimens measuring 6.4 mm in diameter and 
0.5 mm in thickness tested using an MTS Model 859 testing machine 
by indentation with a custom-built, hemispherical head punch at three 
different punch displacement rates (0.5mm/min, l.Omm/min, and 
5.0 mm/min). The tests were performed at room temperature 
(20°C±2°C), according to ASTM F21 83-02 [7]. Again, a total of 
60 specimens were tested (5 specimens per material x 4 materials 3 
rates). 

The experimental data obtained from the small punch tests were in 
terms of force and displacement of the punch head. These data can be 
converted to stress and strain by simulating the small punch tests using 
the finite element (FE) method. FE simulations were performed using 
an axisymmetric representation of the experimental setup (Fig. 1). The 
test specimen was modeled with 360 quadratic axisymmetric elements, 
and the punch and fixture were assumed to be rigid. The FE simula- 
tions were performed using a commercial FE package (ABAQUS, 
HKS Inc., RI). 

An accurate constitutive model capable of predicting the material 
behavior in large-strain multiaxial deformation states is required to 
determine the stress and strain distributions in the test specimens. 
We used the HM [9-11], which has been calibrated and validated 
for the UHMWPE materials in this study. The material parameters 
that were used and a brief summary of the HM are presented in 
Appendix A. 

To enable failure predictions of the small punch testing, a total of 
12 FE simulations were performed (4 materials x 3 rates). In each sim- 
ulation the punch was displaced to the maximum displacement value 
that was experimentally observed in the five repeated experiments on 



that specific material and applied displacement rate. The maximum 
stress and strain values at the experimentally observed failure states 
for each of the five repeated tests were then extracted from the FE 
simulations by finding the location having the maximum stress and 
strain and extracting the stress and strain tensors at that location. 

Failure modeling 

The failure parameters for each material (e.g., the Mises stress at 
failure) were obtained from the average of the experimentally observed 
values from all uniaxial and small punch tests. The quality of the dif- 
ferent failure models was evaluated by comparing the predicted failure 
values for the two loading modes at each of the three displacement 
rates. The quality of the predictions was quantified with the coefficient 
of determination, r 2 , between the experimental failure values and the 
predicted failure condition for each material [9]. The predicted failure 
condition was assumed to be independent of loading rate. 

The predictive capability of each model was also graphically repre- 
sented using “box and whiskers” plots. In these plots, the center 50% 
percentiles of the data are represented using a box. The median is rep- 
resented with a horizontal line inside the box. The 90% and 10% per- 
centiles are represented with vertical bars extending outside the box, 
and any outlying data points outside the 90% and 10% percentiles 
are individually represented by small circles. Thus, the relative ability 
of the eight failure models to predict failure was examined by compar- 
ing the predicted failure value in the different experimental tests for 
each of the materials and by assessing if the predicted failure values 
for each material tended to have a constant value at the observed fail- 
ure instances. 

To evaluate the different failure models, the stress and strain distri- 
butions in the test specimens must be known as a function of the ap- 
plied deformation. For the uniaxial tension tests, the distributions 
were directly obtained from the experimentally determined force- 
displacement curves. For the small punch tests, we obtained the 
stress-strain response by simulating the experimental tests using the 
FE simulations based on the Hybrid Model (HM) [5,8,9] to represent 
the material behavior. 



Results 

The capabilities of the Hybrid Model to predict the 
behavior of UHMWPE can be seen from a direct 
comparison between experimental data for the highly 
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erosslinked material that was annealed at 110°C and 
predictions from the HM (Fig. 2). The uniaxial tension 
and force-displacement responses show that failure 
points occurred at different applied deformations, but 
the experimental force-displacement responses were all 
very repeatable (Fig. 3). As discussed by Bergstrom 
et al. [9-1 1], the HM is an accurate constitutive model 
for predicting the behavior of the different types of 
UHMWPE investigated in this work (Fig. 3). The fail- 
ure strengths of the two highly erosslinked materials in 
uniaxial tension were lower than for the two conven- 
tional materials, but the highly erosslinked materials 
failed at higher forces than the conventional materials 
when tested in small punch (equibiaxial) loading. This 
behavior is caused by different evolutions of the micro- 
structural deformations when tested in uniaxial and 
multiaxial tests [11]. 

The results for the maximum principal stress, Mises 
stress, and Tresca stress models show that these three 
failure models behaved similarly (Table 1) and had 
almost identical failure conditions (Fig. 4a-c). They also 
showed good predictive agreement between the uniaxial 
and the multiaxial failure data for all materials except 
the highly erosslinked material that was heat treated at 
150 °C, For this latter material, the effective stress value 
was significantly higher in the small punch tests than in 




Fig. 2. (a) Uniaxial tension to failure data, (b) small punch loading to 
failure. The experimental tests were performed on UHMWPE 
prepared by four different thermal and radiation treatments. 





Fig. 3. Comparison between the experimental data and Hybrid Model 
predictions: (a) uniaxial tension to failure, (b) small punch testing to 
failure. The experiments were performed on highly erosslinked 
UHMWPE annealed at 110°C. The failure points are indicated by 
the symbol □. 

the uniaxial tests. A large scatter occurred in the effective 
stress values at failure for these three failure models 
(Fig. 4a-c). The variation in stress at failure for the unir- 
radiated UHMWPE tested in uniaxial tension to failure, 
for example, varied between 150 and 320 MPa. 

The predicted hydrostatic stress in the specimen at 
the time of failure was significantly higher in the small 
punch tests than in the uniaxial tests (Fig. 4d). Hence, 
the hydrostatic stress model was not useful as a general 
failure predictor of UHMWPE. The Coulomb model 
was both qualitatively and quantitatively similar to the 
Mises, Tresca, and maximum principal stress models 
(Fig. 4e). For all materials, the maximum principal 
strain at failure was larger in uniaxial tension than in 
the small punch experiments (Fig. 4f). This is reasonable 
since the strain state at the tip of the small punch 
specimen at the time of failure is close to biaxial and typ- 
ically lower than in uniaxial tension. The experimentally 
determined Mises strain at failure was substantially dif- 
ferent in uniaxial tension and small punch loading, 
clearly indicating that the Mises strain was not an accu- 
rate general failure predictor (Fig. 4g). 

In contrast to the other failure models, using the 
chain stretch failure model provided good agreement 
between the predicted failure in uniaxial tension and in 
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Table 1 

Experimentally determined failure conditions for the eight failure models and the four types of UHMWPE 



Failure condition 


Unirradiated 


30 kGy (y-N 2 ) 


100 kGy (110 °C) 


100 kGy (150 °C) 


Max principal stress (MPa) 


246.2 (39.7) 


204.1 (25.6) 


172.1 (25.4) 


158.2 (33.8) 


Mises stress (MPa) 


242.1 (39.1) 


200.6 (25.6) 


170.3 (24.4) 


156.6 (31.3) 


Tresca stress (MPa) 


246.9 (39.9) 


204.4 (25.7) 


172.6 (25.6) 


158.6 (34.0) 


Hydrostatic stress (MPa) 


120.8 (44.8) 


99.2 (33.5) 


85.9 (35.4) 


81.1 (39.6) 


Coulomb stress (MPa) 


271.8 (44.3) 


224.8 (28.2) 


190.1 (28.6) 


174.4 (37.4) 


Max principal strain 


1.48 (0.18) 


1.37 (0.19) 


1.12(0.14) 


1.13 (0.09) 


Mises strain 


1.61 (0.80) 


1.46 (0.70) 


1.12(0.60) 


1.19 (0.59) 


Chain stretch 


2.90 (0.22) 


2.60 (0.21) 


2.04 (0.08) 


2.01 (0.08) 



small punch loading, for all four materials (Fig. 4h). 
Further, the variability in the critical chain stretch be- 
tween different tests performed at different loading rates 
was small relative to the other failure models. The value 
of the limiting chain stretch (A^ ck ), a material parameter 
in the HM that specifies how much the molecules can 
stretch before becoming fully extended, reflected that 
failure occurred when the molecular chain stretch 
reached a certain fraction (60-80%) of the fully extended 
chain stretch (Fig. 4h, Table 2). 

A direct comparison between the quantitative per- 
formances of the different failure models, by examina- 
tion of the revalues between all experimental failure 
data and the critical failure values for the four materials, 
demonstrated that the predictions from the chain stretch 
model had lower variability between different tests and 
better overall agreement with both uniaxial and multiax- 
ial experimental failure data than the other failure mod- 
els (Fig. 5). The chain stretch model described 85% of 
the variability in the failure data for a broad range of 
UHMWPE materials, loading conditions, and loading 
rates. Of the other models, the Mises stress, the maxi- 
mum principal stress, Tresca stress, Coulomb stress, 
and maximum principal strain were about equal, each 
with an r 2 - value of about 0.55. The hydrostatic stress 
and the Mises strain models were clearly the poorest per- 
formers, both with revalues less than 0.2. 



Discussion 

Our study demonstrates that uniaxial and multiaxial 
failure of conventional and highly crosslinked 
UHMWPE under monotonic loading can be accurately 
predicted using the chain stretch model. The model is 
interesting not only because it provides the most accu- 
rate predictions of the models tested, but also due to 
its relationship to the stretching of the macromolecules 
in the material and with the maximum chain locking 
stretch parameter in the Hybrid Model. The correlation 
between the experimental chain stretch at failure and the 
chain locking parameter suggests that if the HM is cal- 
ibrated for a given UHMWPE material, then a good 
multiaxial failure condition can likely be obtained by 



taking the chain stretch at failure to be about 60-80% 
of the chain locking stretch parameter. 

In practical terms, the chain stretch at failure for 
UHMWPE materials can be directly measured during 
standard uniaxial tensile tests, since it is related to the 
uniaxial strain at failure (e): 




Consequently, the chain stretch failure criterion outlined 
in the present study can readily be calculated for con- 
temporary UHMWPE materials based on available 
experimental data. 

Some of the other failure models, e.g. the maximum 
principal stress and the Mises stress models, provide rea- 
sonable failure predictions of the uniaxial and small 
punch data. The scatter in the experimental failure 
criteria from these models, however, makes them less 
useful as a failure predictor under multiaxial loading 
conditions. 

Our results also suggest that the failure conditions for 
conventional and highly crosslinked UHMWPE depend 
on radiation dose and heat treatment (Table 1). For 
example, the failure of the unirradiated and the sterilized 
materials are both nicely captured by either a stress- 
based model (maximum principal stress, Mises stress, 
or Tresca stress) or by the chain stretch model. How- 
ever, for the two highly crosslinked materials (particu- 
larly the remelted material 100 kGy, 150 °C), the 
stress-based failure models tend to underestimate the 
strength in the small punch experiments and overesti- 
mate the strength in the uniaxial experiments. The only 
model capable of accurately predicting the behavior of 
these materials is the chain stretch model. 

It is also interesting to consider the molecular weight 
between crosslinks and how it may relate to the maxi- 
mum chain stretch. The concept of a molecular weight 
between crosslinks is derived from Flory’s [14] statistical 
theories for a swollen crosslinked polymer network in 
equilibrium with a good solvent. In contrast, the maxi- 
mum chain stretch in our failure theory reflects the 
amount of local deformation that can be accommodated 
by UHMWPE, regardless of whether or not it is cross- 
linked. For virgin UHMWPE, the molecular weight 
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Fig. 4. Predicted failure conditions in uniaxial tension and small punch loading of the four different types of UHMWPE examined: (a) maximum 
principal stress; (b) Mises stress; (c) Tresca stress; (d) hydrostatic stress; (e) Coulomb stress; (f) max principal strain; (g) Mises strain; (h) chain stretch. 



between crosslinks is theoretically infinite; thus, there is 
no relationship between molecular weight between cros- 
slinks and maximum chain stretch for the unirradiated 



UHMWPE examined in our study. With respect to 
our materials, the molecular weights between crosslinks 
spanned from infinity (for unirradiated UHMWPE) to 
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Table 2 



Comparison between the experimentally observed chain stretch at 
failure and the chain locking stretch parameter 





Unirradiated 


30 kGy 
(y-N 2 ) 


100 kGy 
(1I0°C) 


100 kGy 
(150 °C) 


Chain stretch 
at failure 


2.90 


2.60 


2.04 


2.01 


Chain locking 
stretch parameter 


4.44 


4.40 


2.80 


2.52 




Fig. 5. Average revalue for all four materials for the different failure 
models. 



an estimated value of 6000 g/mol at a dose of 100 kGy 
[26]. For this reason, we do not judge the molecular 
weight between crosslinks to be an attractive state vari- 
able for predicting the failure behavior of conventional 
and highly crosslinked UHMWPE. 

The material response of UHMWPE is rate-depend- 
ent [9]. A direct analysis of variance of the experimental 
data at different loading rates showed that the only sta- 
tistically significant (p < 0.05) rate-effects in the failure 
quantities were that the punch force at failure increased 
with increasing punch rate and the punch displacement 
at failure decreased with increasing punch rate for 
the unirradiated material. The other failure quantities, 
e.g. the stress and strain response in uniaxial tension, 
showed no significant rate-dependence [22]. Overall, 
these findings support our assumption that failure crite- 
ria for UHMWPE are effectively rate independent. 

The high stress and strain values at failure (Table 1) 
are indicative of the intrinsic strength of the material 
and can be used to predict material failure on the micro- 
scale at the tip of microcracks or the local damage dur- 
ing wear. A complicating factor when studying failure of 
polymeric components is that failure of seemingly iden- 
tical components typically occurs at different load levels. 
The apparent distribution in failure strengths is caused 
by variations in the internal state of the material in dif- 
ferent components. The differences can, for example, be 
due to localized differences in molecular weight, molec- 
ular weight distribution, crosslinking density, degree of 



crystallinity, morphology, or defect population. The 
failure on a macroscopic scale therefore often appears 
stochastic. The results presented here for monotonic 
loading lay the foundation for the future development 
of more advanced stochastic damage models that can 
accurately predict failure under multiaxial, cyclic load- 
ing conditions. 
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Appendix A 



This appendix briefly summarizes the governing 
equations of the Hybrid Model (HM). For a more in 
depth discussion, see the recent publications by Berg- 
strom et al. [9,10]. 

The kinematic framework used in the HM is based on 
a conceptual homogenization of the semicrystalline and 
amorphous domains [3,4,16] and a decomposition of the 
deformation gradient into elastic and viscoplastic com- 
ponents: F = F 6 F p . The Cauchy stress acting on the sys- 
tem can be obtained from the linear elastic relationship: 

T = (2^ e E e + xtr[E e ]l), (A.l) 

where E e = ln[V ®] is the logarithmic true strain, f = 
det[F e ] and /t e , X e are Lame’s constants. For convenience, 
the Lame constants can be replaced by the equivalent 
Young’s modulus ( E) and Poisson’s ratio (v): 

EJ( 2(1 + v e )) and X c = E e vJ(( 1 + v e )(l - 2v e )). 

The details of the stress expressions are: 

Ta = T— 1 — [Tgc ha in(F p ; fi A , #*, ka) + f A T,,(F; /i A )], 

1 + <7 a 

(A.2) 

T — — L 1 (^ P /^lock) ffcy[R P 1 /A 



T/j = Ha 



j tp* 

/p*B p * --+-1 - (IF”) 2 



(A.4) 



where T A is a tensor-valued function of the viscoplastic 
deformation gradient F p and the material parameters 
{^ A ,2 A ck , K A ,q A }, and where p A is the shear modulus, 
2 A ck is the locking stretch, k a is the bulk modulus, q A 
is a material parameter specifying the relative magni- 
tudes of T 8chain and T/ 7 , and B p * is the left Cauchy- 
Green deformation tensor. 

The stress driving the viscoplastic flow of the back- 
stress network is obtained from the same hyperelastic 





374 



J.S. Bergstrom et al / Journal of Orthopaedic Research 23 (2005) 367-375 



Table 3 

Hybrid Model material parameters for each of the four UHMWPE materials used in the tension and small punch simulations 



Material parameter 


Unirradiated 


30 kGy (y-N 2 ) 


100 kGy (110°C) 


100 kGy (150 °C) 


E e (MPa) 


1574 


2020 


2009 


1270 


V e 


0.46 


0.46 


0.46 


0.46 


Ha (MPa) 


7.73 


8.22 


10.15 


8.14 


; lock 


4.44 


4.40 


2.80 


2.52 


k a (MPa) 


2000 


2000 


2000 


2000 




0.20 


0.20 


0.20 


0.20 


*Bi 


40.0 


40.0 


40.0 


40.0 




10.0 


10.0 


10.0 


10.0 


Pb 


27.0 


27.0 


27.0 


27.0 


rg ase (MPa) 


26.1 


25.0 


26.2 


20.7 


m B 


9.5 


9.5 


9.5 


9.5 


tJ??* (MPa) 


8.0 


8.0 


8.5 


8.5 


me 


3.3 


3.3 


3.3 


3.3 



representation that was used to calculate the backstress: 
T b = • Ta(F^), where s B is a dimensionless material 

parameter specifying the relative stiffness of the back- 
stress network. The parameter s B evolves with imposed 
plastic deformation to capture the distributed yielding: 
s’b = * (‘S’b — SBf ) • y v c , where p B is a material parameter 

specifying the transition rate of the distributed yielding 
event, s B j is the final value of s B reached at fully devel- 
oped plastic flow, and y B is the rate of viscoplastic flow 
of the time-dependent part of network B: 

= 1 (A- 5 ) 

where z B = ||dev[T B ]|| F , tg ase and m B are material param- 
eters, and y 0 is a constant coefficient with a value of 1/s. 
The velocity gradient of the viscoelastic flow of the 
backstress network is given by 



f B F e B - 



i dev[T B ] 

tb 



K 



(A.6) 



The yielding and plastic flow of the material is captured 
using 

L p = 7c 



ReT dev[T c ] 



R e 



(A. 7) 



where L p = F P F p -\ T C =T- [F e (T A + T„)F eT ]/./ e is the 
stress acting on the relaxed configuration convected to 
the current configuration, t c = ||dev[T c ]||F is the effec- 
tive shear stress (calculated using the Frobenius norm) 
driving the viscoplastic flow, y c = y o * ( T cAc ase ) mc is 
the magnitude of the viscoplastic flow, r £ ase and m c 
are material parameters, and y 0 a constant coefficient 
with a value of 1/s. 

In total, the augmented HM contains 13 material 
parameters: 2 small strain elastic constants (E e , v e ); 4 
hyperelastic constants for the back stress network 
(p A , 2 l ^ ck , /c A , #a)j 5 flow constants of the backstress net- 
work (s B i,SBf,/^ Tg ase , w B ); and 2 yield and viscoplastic 
flow parameters (r c ase ,m c ). Of these 13 material para- 



meters, four (elastic modulus (£); yield strength (i£ ase ); 
the effective stiffness after yield (ju A ); and the limiting 
chain stretch (A^ ck ), which controls the large strain 
behavior) have been shown to vary between UHMWPE 
formulations, while the remaining nine are constant. 
The material constants used in the analysis are listed 
in Table 3. The HM as presented here has, to date, been 
validated for UHMWPE formulations within the crys- 
tallinity range of 46-61% and for both uncrosslinked 
and crosslinked (100 kGy y-N 2 ) UHMWPE. 
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